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Background The spatial distribution of charge and magnetization within the proton is encoded in the elastic form factors. 
These have been precisely measured in elastic electron scattering, and the combination of proton and neutron form factors 
allows for the separation of the up- and down-quark contributions. 

Purpose In this work, we extract the proton and neutron form factors from world’s data with an emphasis on precise new 
data covering the low-momentum region, which is sensitive to the large-scale structure of the nucleon. From these, we 
separate the up- and down-quark contributions to the proton form factors. 

Method We combine cross section and polarization measurements of elastic electron-proton scattering to separate the proton 
form factors and two-photon exchange (TPE) contributions. We combine the proton form factors with parameterization 
of the neutron form factor data and uncertainties to separate the up- and down-quark contributions to the proton’s 
charge and magnetic form factors. 

Results The extracted TPE corrections are compared to previous phenomenological extractions, TPE calculations, and direct 
measurements from the comparison of electron and positron scattering. The flavor-separated form factors are extracted 
and compared to models of the nucleon structure. 

Conclusions With the inclusion of the precise new data, the extracted TPE contributions show a clear change of sign at low 
Q^, necessary to explain the high-Q^ form factor discrepancy while being consistent with the known ^ 0 limit. We 
find that the new Mainz data yield a significantly different result for the proton magnetic form factor and its flavor- 
separated contributions. We also observe that the RMS radius of both the up- and down-quark distributions are smaller 
than the RMS charge radius of the proton. 

PACS numbers: 25.30.Bf, 13.40.Gp, 14.20.Dh 


I. INTRODUCTION 

The electromagnetic form factors of the proton and 
neutron, and G^’"\ are fundamental quantities 

which provide information on the spatial distributions of 
charge and magnetization in nucleons. The form factors 
are measured using electron scattering where the incident 
electron scatters from a nucleon target through the ex¬ 
change of a virtual photon which serves as the sole medi¬ 
ator of the electron-nucleon electromagnetic interaction. 
By increasing Q^, the four-momentum squared of the vir¬ 
tual photon, the virtual photon becomes more sensitive 
to the small scale internal structure of the nucleon. 

In electron scattering there are primarily two meth¬ 
ods used to extract the proton form factors. The first is 
the Rosenbluth or Longitudinal-Transverse (LT) separa¬ 
tion method [I[, which uses measurements of the unpo¬ 
larized cross section, and the second is the polarization 
transfer or polarized target (PT) method [2|, which re¬ 
quires measurement of the spin-dependent cross section. 
A significant difference is observed between LT and PT 
extractions of the proton form factors ii , which is cur¬ 
rently believed to be the results of larger-than-expected 
two-photon exchange (TPE) contributions [^, @ . 

In this paper, we extract the proton form factors from 
a combined analysis of LT and PT measurements, ac¬ 
counting for TPE contributions following the approach 
of Ref. [Tj , but with an emphasis on recent low-Q^ data. 


The combined analysis allows us to extract the TPE con¬ 
tribution to the e-p elastic cross section, which we com¬ 
pare to other phenomenological extractions and to di¬ 
rect calculations of TPE effects meant to be valid in the 
lower regime, as well as to recent direct measurements 
of TPE contributions from the comparison of electron- 
proton and positron-proton elastic scattering [1,0 . 

With the inclusion of neutron form factor measure¬ 
ments, we separate the nucleon form factors into their up- 
quark and down-quark contributions ©[HI- New cross 
section and polarization measurements at low allow 
for a more detailed examination of the flavor-separated 
form factors in the region sensitive to the large-distance 
behavior of the proton charge and magnetization distri¬ 
butions. We compare these results to previous extrac¬ 
tions, as well as to models of the up- and down-quark 
contributions to the proton form factors. 


II. OVERVIEW OF THE METHODS 

A. Rosenbluth Separation Method 

In the Rosenbluth separation method, the reduced 
cross section an for elastic e-p scattering in the Born 
or one-photon exchange (OPE) approximation is: 

a« = (G^(g2))2 + £(G|(Q2))2^ (1) 
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where t = /AMMp is the mass of the proton, and e 
is the virtual photon longitudinal polarization parameter, 
defined as = [l + 2(1 + t) tan^(^)], where 0e is the 
scattering angle of the electron. For a fixed value, the 
reduced cross section is measured at several £ points, 
and a linear fit of (Ta to e gives as the intercept 

and (G^)^/r as the slope. 

By assuming isospin and charge symmetry and neglect¬ 
ing strange quarks contribution, the nucleon form factors 
can be expressed in terms of the up- and down-quark con¬ 
tributions [Sill 


G n _ ^ 

E,M — '^^E,M 
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E,M^ 


( 2 ) 


where in this convention G^ represents the contribu¬ 
tion from the up-quark distribution in the proton and the 
down-quark distribution in the neutron. 

Solving for G^ ^ and G% from Eq. Q above, we 
get the following expression for the up- and down-quark 
contribution to the proton form factors 
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In the limit = 0, this yields G^ = 2, G% = 1, G^ = 
(2|Ip + ^J.n) = 3.67/rAr, and G^ = (/rp-|-2^„) = -I.OS/xat, 
where /xat is the nuclear magneton. 


B. Recoil Polarization Measurements 


In the recoil polarization method, a beam of longitu¬ 
dinally polarized electrons scatters elastically from un¬ 
polarized proton target. The electrons transfer their po¬ 
larization to the unpolarized protons. By simultaneously 
measuring the transverse, Pj, and longitudinal, Pj, po¬ 
larization components of the recoil protom one can de¬ 
termine the ratio MpG^/G^ in the OPE @,[il[3: 


R = 



Pt {E + E') 
Pi 2Mp 


tan(y). 


( 4 ) 


where E and E' are the initial and final energy of the in¬ 
cident electron, respectively. The ratio can be extracted 
in a similar fashion using polarized beams and targets 
by measuring the asymmetry for two different spin direc¬ 
tions H, 13- 

The two methods yield strikingly different results, with 
values of ^pG^/G^ differing almost by a factor of three 
at high Q^. In the LT separation method, the ratio shows 
approximate form factor scaling, /ipG^/G^ Ri 1, albeit 
with large uncertainties at high values. The recoil 
polarization method yields a ratio that decreases roughly 
linearly with increasing Q^, with some hint of flattening 
out above 5 (GeV/c)^. 


C. Two-Photon Exchange Contributions 


To reconcile these measurements, several studies sug¬ 
gested that missing higher order radiative corrections to 
the electron-proton elastic scattering cross sections, in 
particular two-photon exch ang e (TPE) diagrams, may 
explain the discrepancy |16l - [la | . The role of TPE effects 
was studied extensively both theoretically [Tol - I^ and 
phenomenologically [la [2§-[^. Most calculations sug¬ 
gested that the TPE corrections are relatively small, but 
have a significant angular dependence which mimics the 
effect of a larger value of G^. Detailed reviews of the 
role of the TPE effect in electron-proton scattering can 
be found in H, @ . 

Experimentally, several measurements were performed 
to verify the discrepancy [s^, and to try and measure 
or constrain TPE contributions. Precise examinations of 
the e dependence of or [H, |2^ found no deviation 
from the linear behavior predicted in the OPE approx¬ 
imation. Another measurement was performed to look 
for TPE effects by extracting M^G^/G^ at fixed as 
a function of scattering angle (35l ]. In the Born approx¬ 
imation, the result should be independent of scattering 
angle, and no deviation from the OPE prediction was 
observed. 

Based on the observations above, it is possible to try 
and extract the TPE contributions based on the observed 
discrepancy between the LT and PT results. Assuming 
that the TPE contributions are linear in e and that the 
PT results do not depend on e, and knowing that the 
TPE contribution must vanish in the forward limit (e —> 
1) it is possible to extract the TPE contribution 

to the unpolarized cross section in a combined analysis 
of LT and PT data 0, [H Where polarization 

data exist as a function of £, it is possible to attempt to 
extract the TPE amplitudes with fewer assumptions [s^, 
1^ . though with relatively large uncertainties. 

The most direct technique for measuring TPE is the 
comparison of electron-proton and positron-proton scat¬ 
tering. The leading TPE contribution comes from the 
interference of the OPE and TPE amplitudes, and so 
has the opposite sign for positron and electron scatter¬ 
ing. The only other first-order radiative correction which 
depends on the lepton sign is the interference between di¬ 
agrams with Bremsstrahlung from the electron and pro¬ 
ton, and this contribution is generally small. Thus, after 
correcting the measured ratio for the Bremsstrahlung in¬ 
terference term, the comparison of positron and electron 
scattering allows for the most direct measurement of TPE 
contributions. 

The ratio R^+e- {e, Q^) is defined as 


praw 

P-e+e- 




a{e'^p —>■ e^p) 
a(e~p e~p) 


( 5 ) 


After correcting for the electron-proton Bremsstrahlung 
interference term and the conventional charge- 
independent radiative corrections, the cross section 
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ratio Re+e- reduces to: 

i?e+e- (e, Q") = ~ 1 - 2^27, (6) 

where 62 -) is the fractional TPE correction for electron- 
proton scattering. Until recently, there was only limited 
evidence for any non-zero TPE contribution from such 
comparisons , as data were limited to low or large 
e, where the TPE contributions appear to be small. In 
addition, the details of the radiative corrections applied 
to these earlier measurements are not always available, 
and it is not clear if the charge-even corrections were 
applied in all cases. New measurements have found 
more significant indications of TPE contributions at low 
e and moderate , which are consistent with hadronic 
TPE calculations (^ . 


D. Extraction of the Form Factors and TPE 

To extract the proton form factors, we assume TPE 
contributions to the polarization data are negligible and 
account for the contribution to (T/j by adding the real 
function F{e, to the Born reduced cross section: 

o’ij = (Gm) [1 +j] + (7) 

T fip 

where R = is the recoil polarization ratio. 

We recently extracted the proton form factors and the 
TPE correction a{Q‘^) for large values B, O using 
the world data on gr from Refs. [H, [s^, 1^^. The 
form factors were extracted based on the parametriza- 
tion from Borisyuk and Kobushkin (BK parametriza- 
tion) [s^, where gr is expressed in the following form 

^r={G\j) [1 -I--I-2a(Q^)(l - e)(G^j) . (8) 

This form accounts for the experimental and theoretical 
constraints presented in Sec. Ill Cl The value for the ratio 
R = fJ'pG^/G^ is taken from a parameterization of the 
recoil polarization data, for which we assume there is not 
TPE contribution. In the analysis of Refs. Bini , the 
linear parameterization of R and its uncertainties from 
Ref. |45l | was used. An additional uncertainty was ap¬ 
plied in the analysis of Ref. m to provide more realistic 
uncertainties at lower but this was not included in 
Ref. 0. 

In this work, we extend the analysis of Ref. [ll| to lower 
values by including data from new, high-precision 
cross section measurements, and by improving the pa¬ 
rameterization of R and its uncertainty from polarization 
measurements at both low and high values. Taking 
the recoil-polarization and polarized-target results from 
Refs. we find that the form factor ratio is well 

parameterized as: 


R = 


I -b 0.1430Q2 _ 0.0086Q4 + 0.0072g6’ 


(9) 


with an absolute uncertainty in the fit given by: 

— (0.006)^ + (0.0151n(l + Q^j)^, with in 
(GeV/c)2. 

An argument could be made that the uncertainty 
might be below 0.6% for the low values, especially 
as we approach = 0 where the ratio is known. How¬ 
ever, TPE calculations on the recoil-polarization ratio 
show effects that can be up to 0.5% [ 2 ^ at low Q^. 
Because we neglect TPE corrections to the polarization 
measurements, the uncertainties should be large enough 
to account for this. Figured] shows our new global fit to 
the world data on the recoil-polarization ratio /j,pG^/G^ 
along with the fit’s uncertainty bands. 



FIG. 1: (Color online) The world data on the recoil- 
polarization ratio RpG^/G^ from Refs, along with 

our new fit (black solid line) and its uncertainty bands (black 
dashed lines) and previous fits [3, 1^ . 1^ . 


III. RESULTS 


A. Form Factors and TPE Contributions 


We extract the form factors and the TPE contributions 
using Eq. dSj), following the procedure of Ref. [HI- We 
obtain the value for R and its uncertainty from the new 
parameterization presented in Eq. ®. The cross section 
data included are those used in Ref . a long with the 
addition of low-Q^ data from Refs. 

Figure d] shows the extracted proton form factors 
from this work; these form factors and the TPE con¬ 
tribution are included in the online Supplemental Mate¬ 
rial jl^. We also show the form factors as extracted af¬ 
ter applyin g ha dronic TPE calculations, labeled “AMT- 
Hadronic” 611 and “VAMZ” Ref. [^. In addition. 
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FIG. 2: (Color online) G^/Gu (top) and G\^/{^.-pGd) (bot¬ 
tom) as obtained using the BK parametrization from the 
data of Refs, [s^, Is^ . In addition, we compare 

our results to extractions from several previous TPE calcu¬ 
lations and phenomenological fits: AMT [^, VAMZ f^ . 
Kelly [^, Lomon-GK [^, ABGG [^, Bernauer [^, Ar¬ 
rington Vi-y and Puckett (the fit labeled “new” in 

Ref. [13]). 


we show fits to the form factors from previous phe¬ 
nomenological analyses: “Kelly” “Lomon-GK” 
“ABGG” [si], “Bernauer” “Arrington-y 27 ” il, 

and “Puckett” [13 . 

While extractions using calculated TPE corrections 
and phenomenological-based fits are in reasonably good 
agreement, extractions which combine cross section and 
polarization results but do not allow for an explicit TPE 
correction, in particular the Kelly and Lomon extrac¬ 
tions, tend to have larger differences. The “Arrington- 
y 27 ” fit is very different at large because it is the only 
analysis that allows for TPE contributions to the recoil 


polarization data, though the extraction of these terms 
is extremely model dependent. 

At low , our extraction of is significantly above 
most previous fits. This reflects the discrepancy between 
the recent Mainz data which yields values of which 
are s^ematically 2-5% larger than previous world’s 
data [^. At low this corresponds to only a small dif¬ 
ference in the cross section at large scattering angle, but 
for the larger values of the Mainz experiment, this cor¬ 

responds to a significant difference in the measured cross 
sections. Note that except for the Bernauer result, most 
of the previous phenomenological extractions of the form 
factors and TPE contributions were focused on large 
values, and so did not always worry about how well the 
parameterizations of R reproduced low data. 

Figure [3] shows the TPE term a(Q^) extracted from 
this work. We also show parameterizations of a{Q^) 
from the TPE hadronic calculations of Ref. (BMT) 
and Ref. (AMT), which adds and additional phe¬ 
nomenological TPE contribution at higher and from 
previous phenomenological extractions [3 [fig ■ At high- 
values, the TPE term is consistent with our previous 
extraction and shows an increase in magnitude with in¬ 
creasing Q^. The previous extraction (3 was not well 
constrained at low Q^, with cross section data limited to 
> 0.39 and a parameterization [13 of the polariza¬ 
tion data that was not well constrained at low Q^. So 
below ss 1 (GeV/c)^, the behavior was driven by the 
fitting function which took a{Q^) proportional to y/Q^. 
In the present work, the inclusion of the Mainz data and 
low polarization provide meaningful uncertainties for 
a{Q^) below = 1 (GeV/c)^ and show a change in 
sign at low as seen in previous low-Q^ TPE calcula¬ 
tions [13,111,163, Hi. 

While the Mainz data show indications of structure 
in a{Q^), the uncertainties for these data are underesti¬ 
mates of the true uncertainties. The quoted uncertainties 
on the individual cross sections do not include correlated 
systematic effects, which are a significant contribution to 
the total uncertainty in their final form factor parame¬ 
terization. We fit to the dependence of the extracted 
TPE contributions, adding an additional uncertainty, 
da = 0.01, to each point to help offset the impact of 
the underestimated uncertainties in our extraction from 
the Mainz data. We obtain a{Q‘^) = 0.016 — 0.030-\/^, 
with in (GeV/c)^. Note that for our parameterization 
of TPE contributions, F(Q‘^,e) / {Gj,^ Y i s a linear func¬ 
tion in e, while the —)> 0 limit [y] [ZO] is more nearly 

linear when taken as a ratio to the reduced cross section, 
{G\jY + {e/ t){G^^)‘^, and so the behavior at small Q^, 
where 1/r becomes large, is very different. This can be 
seen more clearly in examining the e dependence of the 
TPE contributions at very low Q^, as shown in Fig. 2] 
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FIG. 3: (Color online) The TPE term a{Q^) as obtained using 
the BK parametrization from the data of Refs. 

Idl . [ 5 ^ . Also shown are a simple parameterization of our re¬ 
sults and curves re pre senting a{Q^) as determined in previous 
analyses [3, [13,1^. l66l| . The bottom plot is on a smaller verti¬ 
cal scale, excluding points with uncertainties above Sa = 0.03 
for clarity, although they are included in our fit. Note that 
we do not show extractions of a{Q^) for the calculations or 
extraction of Ref. [oH] at very low , as the £ dependence is 
quite different in our parameterization when the cross sections 
is dominated by the charge form factor (see text). 


B. Ratio of Positron-Proton to Electron-Proton 
Elastic Scattering Cross Sections Rg+g- 


From the TPE contributions extracted based on the 
parameterization of Eq. ([5]), the corrected ratio of 
positron-to-electron scattering cross sections is 


Re+e- 


1 — ^27 
1 -l- (^27 


1 - 


4a((5^)(l - e) 

(1 + ^) 


( 10 ) 


Figure |4] shows the ratio Re+e- as a function of e ex¬ 
tracted for a range of values. The solid red curve 
represents the ratio from our new parameterization of 
a((5^), and the others show the results of TPE hadronic 
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FIG. 4: (Color online) The ratio R^+e- {£,Q^) as determined 
from our parameterization of a(Q^) at the value listed in 
the figure. The other curves show previous calculations [6l|] or 
phenomenological extractions [3. [^. T he data points 

are direct measurements of Rg+g- [3,[3,[zII“E3]- For the world 
data, the measurement and value(s) in (GeV/c)^ are given. 


calculations from Refs. l20l l6ll a nd several phenomeno¬ 
logical fits from Refs. I'd. 123. 11 in. The data points 
are world data on the ratio Rg+g- from Refs. Silnl- 
IH- Two recent measurements [1, H have provided sig¬ 
nificantly more precise measurements at 1.0 and 

1.5 (GeV/c)^, providing evidence for non-zero TPE at 
larger values and a change of sign from the exact cal¬ 
culation at = 0, consistent with what we observe. Our 
results are slightly larger than the direct measurements 
at 1 (GeV/c)^ from Ref. i, but otherwise in very good 
quantitative agreement with existing data. 


Note that parameterizations where the TPE contribu- 
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tion is similar to that of Eq. (jS)), i.e. a linear function 
times the low results for R^+e- will have a 

strong non-linear behavior. The TPE contribution rela¬ 
tive to is linear, but (G^)^ dominates the cross 

section at very low except for £ 0, strongly sup¬ 

pressing TPE as a fractional contribution as one moves 
away from e = 0. This could be fixed by modifying the 
functional form, e.g. using a linear function in e times the 
full reduced cross section However, because the re¬ 
sults are similar everywhere except at very low , where 
there is little data, we consider this parameterization suf¬ 
ficient for the present analysis. 


C. Flavor Separation of the Nucleon Form Factors 

Examinations of the flavor-separated form factors of 
the nucleon focusing on high Elm have provided 
several interesting observations. Below we summarize 
some of the main conclusions from these analyses: 

(1) The down-quark contributions to Dirac, Fi, and 
Pauli, F 2 , form factors deviate from the expected l/Q'^ 
scaling with small differences between the de¬ 
pendence in Fi and F 2 for both the up- and down-quark 
contributions. 

(2) The up- and down-quark yield very different con¬ 

tributions [ll| to Gb/Gm and F^jF^. The strong linear 
falloff with in the ratio is not seen in ei¬ 

ther the up- or down-quark contributions, but mainly 
arises due to a cancellation between a weaker depen¬ 
dence for the up-quark and a negative but relatively Q^- 
independent contribution from the down-quark. 

(3) The more recent analysis m shows some differ¬ 
ences from the original work by Cates et al. 0 , referred 
to as “CJRW” throughout the text. The difference are 
associated with the approximations made in the CJRW 
analysis, which neglected TPE effects and included only 
the uncertainty associated with the neutron charge form 
factor. The treatment of TPE contributions yields small 
but clear differences, mainly at lower values, up to 
Ri 1.5 (GeV/c)^, while the addition of uncertainties asso¬ 
ciated with TPE and all of the form factors yields some¬ 
what larger uncertainties in most form factors, and pro¬ 
vides an estimate for the uncertainties in the magnetic 
form factors absent in the CJRW analysis. 

Figures [5] and [6] show the flavor-separated form factors 
extracted from this work; the extracted values are in¬ 
cluded in the online Supplemental Material ■ The up¬ 
dated parameterization of the polarization-transfer mea¬ 
surements of fjLpG^^/G\i yields a small difference in the 
data above 5 (GeV/c)^, but it is always a factor of 2-3 
below the quoted uncertainties. Note that because the 
GJRW analysis includes only the uncertainties from G^, 
the flavor-separated extraction of Gm is quoted without 
uncertainty. In our analysis, uncertainties are included 
for all form factors, but we use a parameterization for the 
uncertainty on G'% and G^, taken from Ref. [ll|. Be¬ 


cause of this, cases where the uncertainty is dominated by 
neutron data yield points with large error bars but small 
scatter between points, as only the proton uncertainties 
are independent for each values. 

The results are compared to the original CJRW re¬ 
sults and two recent extractions of the nucleon elas¬ 
tic form factors which apply calculated TPE contribu¬ 
tions [ill, li^- Also shown are calculations based on 
dressed-quarks contribution within the frame work of 
Dyson-Schwinger equation “DSE” from Ref. 8l|, pion- 
cloud relativistic constituent quark model “PC-RCQM” 
from Ref. [s^ , relativistic constituent quark model whose 
hyperfine interaction is derived from Goldstone-boson ex¬ 
change “GBE-RCQM” from Ref. and generalized 

parton distributions with incorporated Regge contribu¬ 
tion that apply diquark models “GPD” from Ref. [s^. 


We start by examining the contributions to G^. The 
down-quark contribution is much smaller than the up- 
quark contribution at all values. Up to Rs 
1 (GeV/c)^, both G’^fGa and G^/Gd increase with 
deviating noticeably from the dipole form at k, 
0.10 (GeV/c)^ and r» 0.40 (GeV/c)^, respectively. 
At higher values, G'^IGo decreases rapidly while 
G%IGd continues to increase, leading to the significant 
falloff in G^^jGo after applying the charge weighting of 
the up- and down-quark contributions. While most of 
the calculations give a reasonable qualitative descrip¬ 
tion of the data, all showing a rise and then fall of 
both G'^/Gd and G'^/Gd with increasing values, the 
AMT-Hadronic and VAMZ parameterizations, as well as 
the GPD model, which also fits to the measured form 
factors, provide the best description of the data. 


This work allows for a more detailed examination of the 
low-Q^ region, which is sensitive to the RMS radius of the 
up- and down-quark distributions. At very low values, 
G'% is consistent with the dipole fit, while G% falls less 
slowly than the dipole, yielding an increase in the ratio 
G%IGd- The proton RMS charge radius as determined 
from electron scattering is te ~ 0.88 fm [sl, [b^. The 
dipole form corresponds to an RMS radius of 0.811 fm, 
suggesting an RMS radius of approximately 0.81 fm for 
the up-quark distribution and below 0.81 fm for the down 
quarks. This gives Vd < < Tp for the charge radii, i.e. 

both the up- and down-quark distributions are more lo¬ 
calized than the overall proton charge distribution, due 
to cancellation of the up- and down-quark contributions 
in the total proton charge distribution. Note that while 
the proton charge radius as extracted from muonic hy¬ 
drogen measurements [H, [s^ yield te ~ 0.84 fm, it is 
more natural to compare scattering results of the proton 
and flavor-separated contributions, although the up- and 
down-quark radii are also smaller than muonic hydrogen 
charge radius. 

For G\,j, the up-quark contribution is again dominant, 
even more so than for the charge form factor. Our re¬ 
sults are in relatively good agreement with the GJRW 
analysis for above 2 (GeV/c)^, though with a small 
offset associated with TPE. For ^ 1 (GeV/c)^, we 
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FIG. 5: (Color online) {G^/Gd) (top) and (G|^d) (bot- 
tom) as a function of from the data of Refs. [33,[3l,[3l- 
13, . Also shown the C JRW extractions [l^ (open trian¬ 

gles), the AMT-Hadronic [hH and VAMZ [62| fits, and the 
values from the GPD (^. DSE [3], PC-RCQM [3], and 
GBE-RCQM [s^ . I^l models. Note that for the proton, the 
up- and down-quark contributions have weighting factors of 
2/3 and -1/3 (Eq. ©). 


FIG. 6: (Color online) {GIi/Gd) (top) and {Gj^ [^^oUhot- 
tom) as a function of from the data of Refs. ill. Slll- 
|3i [Hi • Also shown the CJRW extractions 111] (open trian¬ 
gles), the AMT-Hadronic [H and VAMZ [62l| fits, and the 
values from the GPD [H, DSE [H, PC-RCQM [H, and 
GBE-RCQM [Ul models. Note that for the proton, the up- 
and down-quark contributions have weighting factors of 2/3 
and -1/3 (Eq. ([2])). 


find noticeably larger values for both the up- and down- 
quark contributions to the magnetic form factor. This 
difference comes from the Mainz data which yields val¬ 
ues of which are systematically larger than previous 
world’s data, as seen in Fig. [5] Thus, it is not surpris¬ 
ing that there is an inconsistency between the extrac¬ 
tion from these data and other measurements, in partic¬ 
ular for the up-quark contribution. For both the flavor- 
separated charge and magnetic form factors, the global 
fits [6l|) IH describe the data (excluding the results from 
the Mainz measurement) well. The models with many 
free parameters [H, tend to do a better job than the 
models with few or no parameters adjusted to match the 


form factor data [8l|, [31, although the GPD parameter¬ 
ization provides a significantly worse description of the 
magnetic form factor. 

Focusing on the very low data, we see that G'^ falls 
more slowly than the dipole, indicating a magnetic ra¬ 
dius smaller than the 0.81 fm associated with the dipole 
form. The magnitude of the down quark contribution 
falls more rapidly with indicating a larger radius, al¬ 
though because Gf^ is negative this corresponds to an 
increase in the ratio. The uncertainty on the proton’s 
magnetic radius is sig nificantly larger than for the elec¬ 
tron radius [55[. [3. [^ . so while the magnetization distri¬ 
bution is clearly larger for the down quarks than for the 

























up quarks, opposite of what is observed for the charge 
radius, it is difficult to determine how these compare to 
the overall proton magnetization radius. 


D. Flavor-dependent contribntions to gr 

Using Eqs. m and m, the reduced cross section in the 
Born approximation can be written as 

= + (n) 

which can be separated into terms coming from just the 
up-quark or down-quark contributions and the up-down 
interference term: 

( 12 ) 

Figure [3 shows the total and flavor-separated contri¬ 
butions to the reduced cross section as a function of e 
for a sample of values. The down-quark term, a^') is 
extremely small, often smaller than the TPE corrections. 
The up-quark aG\ and the up-down quark interference 
uG')x(d) dominant contributions, with the up- 

down interference term being comparable in size to the 
total reduced cross section at low values. So while 
down-quark contribution, representing the form factor 
one would obtain if only the down quark distribution was 
present, is almost negligible at all values of £ and Q^, the 
net impact of the down-quarks on the cross section is still 
large, especially at low values. Note that the strange- 
quark contributions have been neglected, and while the 
contributions coming from (G|; will have a negligible 
contribution to the cross section, the up-strange interfer¬ 
ence term will be significantly larger, and could have a 
non-negligible contribution. 

To simplify the examination of the flavor-dependent 
contributions, we examine the contribution of the up, 
down, and interference terms to (G^)^ and (G^)^, as 
shown in Fig. [S] For e = 0, the cross section depends on 
iG%Y which is dominated by the up-quark contribution, 
with a roughly 20% contribution from the interference 
term adding to (G^)^ at all values. 

The breakdown of the charge form factor is more com¬ 
plicated. At low values, the interference term is 
roughly half the size of the up-quark contribution and 
of the opposite sign, making the interference term nearly 
identical in magnitude to the total value of (G^)^. In¬ 
creasing from Ri 0.1 to 1 (GeV/c)^, the decrease in 
G^ is driven by the increase in the magnitude of the 
(negative) up-down interference term, moderated by a 
slight increase in the up-quark contribution. Above 1-2 



e e 



FIG. 7: (Color online) The reduced cross section gr including 
TPE (solid line) and without TPE (dotted black line), along 
with the flavor-separated contributions to the Born cross sec¬ 
tion: gG'> (up-quark), (down-quark), and g^^^G) ('^p_ 
down interference). 


(GeV/c)^, the rapid fall of (G^)^, partially cancelled by 
an increase in both the down-quark and up-down inter¬ 
ference terms, continues the very nearly linear decrease 
observed in the /r^G^/G^ ratios. Thus, the observed 
nearly linear behavior comes from a complicated combi¬ 
nation of the up, down, and interference terms, with none 
of the individual contributions showing such a monoton- 
ically decreasing behavior relative to the dipole form. 

The faster falloff of the down-quark contributions was 
interpreted in Ref. 0 and references therein as an indi¬ 
cation of the possibility of sizable nonzero strange matrix 
elements at large or the importance of diquark de¬ 
grees of freedom. While existing measurements of parity- 
violating elastic scattering yield very small contributions 
from the strange quarks up to ss 1 (GeV/c)^ 
they still leave open the possibility for significant con¬ 
tributions from G|; and GU which cancel in the parity- 
violating observables [oll, , although there are also re¬ 
sults from lattice QGD that the strange-quark contri¬ 
bution is small for the charge and magnetic form fac¬ 
tors (^.1^. 


In the diquark model, the singly occurring down-quark 
in the proton is more likely to be associated with an 
axial-vector diquark than a scalar diquark, and the con¬ 
tribution of the axial-vector diquark yields a more rapid 
falloff of the form factors. The up quarks are generally 
associated with the more tightly bound scalar diqu arks, 
yielding a harder form factors EaSilMlai-Iiol- Re- 
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FIG. 8: (Color online) Flavor-separated contributions to 


cent calculations suggest that pion loop corrections 
play a crucial rule for ^ 1-0 (GeV/c)^. The nucleon 
form factors were expressed as proton quark sector form 
factors, and was found that the down-quark sector of 
the Dirac form factor is much softer than the up-quark 
sector as a consequence of the dominance of scalar di¬ 
quark correlations in the proton wave function. On the 
other hand, the up-quark sector of the Pauli form factor 
is slightly softer than in the down-quark sector suggesting 
that the pion cloud and axial-vector diquark correlations 
dominate the effect of scalar diquark correlations leading 
to a larger down-quark anomalous magnetic moment and 
a form factor in the up-quark sector that is slightly softer 
than in the down-quark sector. 


IV. CONCLUSIONS 

In conclusion, we improved on and extended to lower 
the previous extraction Q of form factors and two- 
photon exchange contributions, as well as the extrac¬ 
tion m of flavor-separated contributions to the proton 
form factors. We used new polarization data to obtain an 
improved parameterization of /ipG^/G^ and its uncer¬ 
tainties, and included new cross section data [H, [H, 
to extend the analysis to lower values. 

The results for G^ are generally in excellent agreement 
with those extracted based on a global analysis including 
calculated TPE hadronic corrections [bll, [^ , as well as 
some previous phenomenological extractions . 

For G^j, the results disagree noticeably with previous 
extractions which applied hadronic TPE corrections and 
as well as some previous phenomenological extractions. 
This is in large part due to the tension between the 
new low-Q^ Mainz data [^, although the different ap¬ 
proaches for applying TPE corrections have some impact 
as well, especially for extractions which neglected TPE. 

Our new low extraction of the TPE contribution 
yields values of a{Q^) at the few percent level. The TPE 
term shows a change of sign at ^ 0.40 (GeV/c)^, 
which was not seen in our previous extractions and 
fit 0, but which is consistent with low TPE calcu¬ 
lations [13, HH in, 113) m • This is the first phenomeno¬ 
logical extractions that directly observes this predicted 
change of sign at very low values. This behavior was 
present in the extraction of Bernauer et al. [q^ ] , but be¬ 
cause the low = 0 limit was fixed to the Feshbach 
correction [t^I) the change of sign is guaranteed by the 
fit procedure. 

We compared our extracted TPE corrections to previ¬ 
ous phenomenological extractions and to world’s data on 
the ratio of positron-proton and electron-proton scatter¬ 
ing cross sections, Re+e- ■ Our extracted TPE contribu¬ 
tions are in generally good agreement with world’s data 
on i?e+e-, including recent measurements which show a 
clear e dependence, consistent with the form factor dis¬ 
crepancy, at values of 1-1.6 (GeV/c)^ ii- 

Inclusion of the new low-Q^ polarization and cross sec¬ 
tion data allows us to examine the slope of the form factor 
at small Q^, which is connected to the flavor contribu¬ 
tions to the RMS charge and magnetization radii of the 
proton. While the low data do not allow for a pre¬ 
cise extraction of these radii, we find that both the up- 
and down-quark distributions are more localized than the 
overall charge distribution, with rd < ru < fp. For the 
magnetization distributions, the up-quark contribution 
has a smaller radius than the down-quark contribution, 
providing another clear indication that the magnetiza¬ 
tion distribution does not simply come from the spatial 
distribution of the quarks. 

We also express the reduced cross section an in terms 
of these flavor-separated from factors to shed light on 
the contributions of the up-quark, down-quark, and up- 
down quarks interference terms to the Born cross section 
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<j. The up-quark contribution dominates at all 
values, while the down-quark contributions repre¬ 
senting the cross section one would observe if scattering 
from only the down quarks, is typically negligible. The 
interference term, is negative and can be size¬ 

able. While it is always smaller than the up-quark con¬ 
tribution, it can be comparable in size to the total cross 
section after the significant cancellation between the in¬ 
terference terms and the up- and down-quark terms. 
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